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Abstract. Models of particle propagation in causal set theory are investigated 
through simulations. For the swerves model the simulations are shown to agree with 
the expected continuum diffusion behaviour. Given the limitations on the simulated 
' causal set size, the agreement is far better than anticipated. 
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1. Introduction 



Causal set theory is a discrete, Lorentz invariant approach to quantum gravity. For 
comprehensive reviews of the field see, for example [H [2j [3]. The phenomenology of 
particles in causal set theory was first investigated by Dowker et al. [4J, and later further 
developed in [5]. Massive particles propagating in a discrete spacetime are expected to 
experience small fluctuations in momentum, an effect termed 'swerves'. Several simple 
classical microscopic models for particle propagation on a causal set were proposed 
in [1] and [5], but it was shown that it was not necessary to choose a specific model to 



understand the phenomenology in the continuum limit. Dowker et al. demonstrated that 
any Lorentz invariant Markovian stochastic process on the massive particle state space 
gives rise to a continuum diffusion equation. In terms of the observers time parameter, 
'cosmic time', the swerves diffusion equation is [5]: 

dp - ~ p ap + kd a Ov^fc ( -V) V (!) 



dt rwy V. 

where p is the probability density for the system to have a certain momentum and 
position at a given time, p is the particle momentum, m is the particle mass, 7 = 
a/ m 2 + p 2 /m is the standard relativistic factor, and g is the metric on the mass shell, 
i.e. the hyperboloid in momentum space defined by p^ = —m 2 . k is the diffusion 
parameter, the one free parameter of the model. 

Although no one microscopic model was chosen to derive the diffusion equation, the 
phenomenological parameter, k, will in reality depend on the properties of the underlying 
model for particles on a causal set. This paper investigates whether the behaviour of 
underlying particle models is well described by the continuum limit swerves diffusion 
equation. As will be shown below, the results are in fact better than expected, the 
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model proposed in [4J is well approximated by the diffusion equation even far from the 
continuum limit. This allows the relationship between the microscopic model parameters 
and the diffusion parameter to be determined. 

First it is necessary to review the basic concepts of causal set theory. A causal set 
is a set C endowed with a binary relation 'precedes', -<, that satisfies: 

(i) transitivity: if x -< y and y -< z then x -< z, Wx, y, z G C; 

(ii) reflexivity: x -< x, Wx G C; 

(iii) acyclicity: if x -< y and y -< x then x = y, Wx, y G C; 

(iv) local finiteness: Wx, z G C the set {y \ x -< y -< z} of elements is finite. 
Some definitions will be useful for the following work. Let C be a causal set. 

(i) A chain is a totally ordered subset of C. 

(ii) A longest chain between two elements x, y G C is a chain whose length is longest 
amongst chains between those endpoints. There may be more than one longest 
chain between two elements. The length of the longest chain between elements 
x, y G C will be denoted d(x, y). 

(iii) A link is an irreducible relation: elements x and y are linked if and only if 
d(x, y) — 1. If two elements x, y G C are linked, it will be denoted x -<* y. 

(iv) A path is a chain consisting of links. 

2. Particle models 

For the purpose of simulations, a causal set can be constructed from a continuum 
manifold by a process called sprinkling. Points are selected at random from a manifold 
via a Poisson process in which the probability measure is equal to the spacetime volume 
measure in fundamental units. The causal order on the points induces the partial order 
on the elements in the causal set. 

Construct a causal set by sprinkling into Minkowski spacetime. A massive particle 
trajectory is taken to be a chain of elements e n in the causal set, i.e. a linearly ordered 
subset of C. The particle trajectory is constructed iteratively. It is assumed that the 
trajectory's past determines its future, but that only a certain amount of the past is 
relevant. 

2.1. Swerves 

First proposed in [3], this model relies on information about the approximating 
spacetime. Suppose the particle is currently located 'on' an element e n , with a four- 
momentum p n . The next element, e n+ i is chosen such that 

• e n+ i is in the causal future of e n and within a propertime 77 of e n , 

• the momentum change \p n +i — p n \ is minimized. 

Tf is the forgetting time of the process. The momentum p n+ \ is defined to be proportional 
to the vector between e n and e n+ i, normalized by the particle mass m. 
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2.2. Intrinsic models 

Two models intrinsic to the causal set, i.e. not relying on continuum information, were 
proposed in [5] to illustrate the range of possibilities available. Rather than a forgetting 
time, these models depend on a forgetting number n/. 

2.2.1. Model 1 Given a partial particle trajectory . . . e„_i, e n the next element e n+ i is 
chosen such that 

• d(e n _i,e n+ i) < 2n f , 

• d(e n , e n+ \) is maximized subject to d(e n , e n+ \) < nj. 

These requirements do not guarantee the existence of a unique e n+ \. There will, however, 
almost surely be finitely many eligible elements and the trajectory can be constructed 
by choosing an element uniformly at random from these. Note that this model is 
slightly different from the first intrinsic model given in [5], where equalities in the above 
conditions were given. Model 1 of [5] does not guarantee the existence of an e n+ i under 
reasonable conditions. 

2.2.2. Model 2 The trajectory is constructed as a path in this model, i.e. d(e n , e n+ i) = 
1 for any e n , e n+ %. Given a partial particle trajectory . . . e n _ n/ , . . . , e n _i, e n the next 
element e n+ \ is chosen such that 

• d(e n ,e n+1 ) = 1, 

• d(e n - nf , e n+ i) + . . . + d(e n - h e n+ i) + d(e n , e n+1 ) is minimized 

Again this minimization does not necessarily yield a unique e n+ i, in which case the 
trajectory is constructed by choosing an element uniformly at random from those 
eligible. Also, if the trajectory has length less than nf the minimization is done over all 
elements available. 

In these models it is assumed that the forgetting parameter t/ or nj is many orders 
of magnitude greater than the discreteness scale. 

3. Numerical results 

Simulations of the particle models given above were developed within the Cactus 
numerical relativity framework [6], making use of the CausalSets arrangement written 
by David Rideout. Although the models are entirely general, due to computational 
limitations simulations were carried out in 1+1 dimensions. Points were sprinkled into 
a region of Minkowski spacetime with a Poisson distribution with a mean number of 
elements N. One additional point was added to each causal set at the origin to give 
a fixed beginning point, ei, for the trajectories. For convenience, the particular frame 
to which the coordinates of the sprinkled points in Minkowski spacetime refer will be 
called the embedding frame. For the swerves model the particle was assumed to be 
initially at rest in the embedding frame, i.e. pi = (m, 0). For the intrinsic models the 
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(a) Model 1, N = 32768, r f = 0.1105, I = 16. (b) Model 2, N = 32768, n f = 20, I = 13. 
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(c) Model 3, N = 32768, n f = 20,1 = 163. 



Figure 1. Example trajectories with causal set size, TV = 32768, forgetting parameter 
Tf/rif and trajectory length I. 



particles were also assumed to be initially close to rest. To construct the first step for 
the intrinsic model 1, the condition d(eo, e<i) < 2nj was neglected and e<i was taken to be 
the first (in time) element that maximized d(ei, C2) < n/. In the case of intrinsic model 
2 an additional point was added to the causal set at x — 0, t — 0.25 and the beginning 
of the trajectory was taken to be a longest chain between the points x = 0, t = 
and x = 0, t = 0.25, to ensure the particle was initially close to rest. An example 
trajectory for each of the models is shown in Figure [TJ Each model gives fluctuations 
in momentum, but the amount of fluctuation is clearly model dependent. Note that 
for a 1+1-dimensional causal set with N = 32768, nj = 20 is roughly equivalent to 
Tf = 0.1105 if we take d p i = y/V/N and n/ = Tf/d p i. 

To investigate how well the diffusion equation approximates the microscopic model, 
and the relationship between the phenomenological parameter and the underlying model 
parameters, we focused on the swerves model. For a given causal set and initial position 
and momentum, the swerves model defines a unique trajectory. To say that the swerves 
model results in diffusion is, in a sense, saying that the exact underlying causal set 
is unknown. To simulate this, many different sprinklings into the same region of 
Minkowski spacetime were generated, and the unique trajectory in each was calculated. 
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For each trajectory, the 'final' position and momentum were determined as the trajectory 
crossed the tf — 0.95 hypersurface. This gave a distribution in position and momentum 
that could be compared to that expected from the 1+1 dimensional swerves diffusion 
equation: 



Before comparing the simulation results and the diffusion equation we must note 
that some trajectories may need to be rejected. If a trajectory is close to the boundary 
of the region of Minkowski spacetime at any point, it will 'bounce' back and distort the 
results. For the results shown here, the models parameters were chosen such that there 
were few such 'invalid' trajectories - any that did occur were removed from the results. 

The microscopic model contains three parameters: the forgetting parameter, Tf, a 
discreteness scale, d p i, and the particle mass, m. Although the discreteness scale does 
not appear in the trajectory algorithm, the trajectory clearly depends on how many 
points have been sprinkled into a given volume. For the purposes of this investigation 
the discreteness length can be defined as d p i = yJV/N, where V is the volume of the 
region of 1+1 dimensional Minkowski spacetime and N is the mean number of causal 
set elements sprinkled. 

The particle mass appears in the swerves model algorithm only to normalize the 
momentum at each step. The trajectory constructed is, in fact, independent of the mass. 
The final momentum distribution from a collection of trajectories does depend on the 
mass, but changing the mass only rescales the momentum. Thus, without running any 
simulations it is possible to determine the dependence of the diffusion parameter, k, on 
m: examining (j2J), if m and p are rescaled by a factor a then k must be rescaled by a 2 , 
therefore k ~ m 2 . 

To determine the dependence of k on Tf, 500 trajectories were evolved for each of 
11 values of Tf between 0.03 and 0.1, with m = 1 and d pt ~ 0.0055, in embedding units. 
For each value of t/ final position and momentum histograms were calculated. The 
swerves diffusion equation was numerically evolved for a range of values of k and a best 
fit value of k was determined for each 77 by minimizing the reduced \ 2 value 



where Oj is the observed frequency for a momentum bin i, Ei is the expected frequency 
(i.e. that given by the evolution of the diffusion equation), and / is the number of degrees 
of freedom (here, / = number of data points - 1). \ 2 is not a good measure of fit if a 
significant proportion of the expected frequencies are less than five (see, e.g. fF]). To 
avoid this problem multiple bins were combined where necessary. A reduced \ 2 of order 
1 is usually considered to indicate a good fit (see, e.g. [8]). The momentum diffusion 
was chosen for the comparison as the position diffusion is driven by the momentum 
diffusion. 






(3) 
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Figure 2. Histograms for 500 trajectories with t/ = 0.04, N = 32768 and best fit 
solutions from the swerves diffusion equation, k — 2.8 x 10~ 9 . The reduced \ 2 value 



for the momentum distribution is x~ 



0.59. 



Example position and momentum histograms and the corresponding best fit 
solutions for Tf = 0.04 are shown in Figure [2j Note that Tf = 0.04 has not been 
chosen for any particular reason, the other values of r/ have equally good fits. It is clear 
from this figure that the swerves diffusion equation is a very good approximation to the 
underlying model even though t/ is not many orders of magnitude greater than d p i and 
the discreteness length is nowhere near the d p i — > continuum limit. 

Plotting Ink vs. lnry, Figure 3(a), for all 11 values of Tf reveals the dependence 
k ~ rf. 

Dimensional analysis now allows the dependence on the final parameter, d p i, to 
be determined, but for completeness it was checked through simulations. For four 
different causal set sizes N = {4096, 8192, 16384, 32768}, 500 trajectories were evolved 
with fixed Tf and m. Again, best fit values of k were determined for each value of N. 
Figure 3(b) clearly shows, despite only four data points, the dependence k ~ d^. Thus, 
k ~ m 2 d p i l /Tf for the swerves model. The constant of proportionality can of course be 
determined: returning to the varying Tf data and working in discreteness units where 



d p i = 1 it is found that k ~ 2m 2 /r|, as shown in Figure 3(c 



This relationship will be crucial if and when there is an independent reason for a 
particular value of Tf. For example, it has been hypothesized that causal set theory 
will contain a 'nonlocality' scale [9], [10]. The forgetting parameter that appears in the 
above particle models may be a measure of this nonlocality scale. Sorkin [9] estimates 
the nonlocality scale to be of the order 10~ 12 cm, or 10 20 in Planck units. If the results 
above are extrapolated to such large values of Tf, a proton with m = 10 -20 would have a 
diffusion parameter of order k ~ 10~ 140 . The diffusion parameter for hydrogen molecules 
(the factor of two mass difference is inconsequential here) has been constrained to be 
k < 10~ 102 |4j. The above estimate is thus not ruled out by existing constraints, but is 
also, unfortunately, too small to be currently tested. 
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(b) In A: vs. lnd p / and best fit line 4.0\nd p l + 13 
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c) k vs. m /t^ and best fit line 1.9 m /rf — 2.5 x 10 
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Figure 3. The relationship between the diffusion parameter and underlying model 
parameters. 
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4. Conclusions 

There is a considerable focus in quantum gravity phenomenology on violations of 
Lorentz invariance despite there being no evidence to date that Lorentz invariance 
is violated. Causal set theory offers a way to investigate Lorentz invariant quantum 
gravity phenomenology. Earlier work led to a diffusion equation describing the behaviour 
of massive particles in a discrete spacetime in the continuum limit. Although some 
underlying models were proposed, no formal connection between the models and the 
continuum behaviour was made and it was expected that the limitations on the size 
of causal sets that can be simulated would prohibit any direct demonstration of the 
diffusion behaviour. The results given here show that microscopic models of particle 
motion do indeed give rise to diffusion and, moreover, it is not necessary to take the 
discreteness scale, d p i, to zero for this behaviour to occur. The models discussed here are 
classical point particle models, and there is no claim that they are the true description 
of particles in causal set theory. This work demonstrates, however, that we can gain 
much useful information about the observable consequences of discreteness through 
simulations, despite computational limitations on causal set size. 
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